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1.  Introduction 


Accurate  measurement  of  in-flight  kinematics  significantly  contributes  to  the  development  of 
experimental  projectiles  and  rockets  and  to  diagnostics  for  existing  munitions  and  weapons 
systems.  Such  measurements  can  in  some  cases  be  made  with  high  speed  photography,  but  this 
technique  is  generally  used  for  only  limited  portions  of  a  projectile  flight  for  reasons  of  expense 
and  practicability.  Also,  the  precision  of  angular  measurements  is  limited  in  this  methodology. 
Another  measurement  technique  used  for  obtaining  angle  of  attack  data  is  yaw  cards,  but  this 
technique  also  is  low  resolution  and  provides  only  a  small  number  of  discrete  data  points  along  a 
trajectory.  Radar  can  provide  position,  velocity,  and  spin  measurements  but  does  not  give 
projectile  orientation  infonnation.  Global  positioning  systems  provide  position  and  velocity 
measurements  but  cannot  provide  orientation  information.  On-board  sensor  systems  fixed  to  a 
projectile  body  combined  with  a  telemetry  system  can  provide  continuous  measurements 
throughout  a  projectile’s  entire  trajectory. 

The  Advanced  Munitions  Concepts  Branch  of  the  U.S.  Army  Research  Laboratory’s  (ARL’s) 
Weapons  and  Materials  Research  Directorate  has  for  many  years  designed,  built,  and  employed 
body-fixed  sensor  systems  in  support  of  military  ordnance  development  and  testing.  Sometimes, 
a  quantity  of  interest  to  a  tester  is  directly  measurable,  and  given  an  accurate  sensor,  only  the 
observation  or  recording  of  sensor  data  is  required  to  meet  test  objectives.  Temperature  is  one 
example  of  such  a  quantity.  In  other  cases,  body-fixed  sensors’  output  needs  to  be  combined 
and/or  otherwise  processed  in  order  to  obtain  measures  of  desired  quantities.  For  example,  a 
post- launch,  time  history  of  both  speed  and  heading  is  required  to  locate  a  flying  body. 
Combining/processing  is  also  often  required  for  correct  interpretation  of  output  from  body-fixed 
sensors  responding  to  stimuli  external  to  the  flight  body.  Vector  magnetometers  measuring 
components  of  the  earth’s  magnetic  field  are  an  example  of  such  sensors. 

In  essence,  there  is  a  twofold  difficulty  in  correctly  interpreting  and  employing  data  from  body- 
fixed  sensors.  First,  sensor  responses  many  times  are  affected  by  stimuli  other  than  those  that  a 
sensor  is  intended  to  quantify,  e.g.,  an  angular  rate  sensor  may  be  affected  by  any  g  forces  to 
which  it  is  subjected.  Second,  sensor  systems  often,  of  necessity,  make  measurements  in  a  body- 
fixed  coordinate  system,  and  the  quantities  whose  values  are  desired  are  best  described  in 
another  coordinate  system.  Although  coordinate  transforms  are  a  staple  of  high  school 
trigonometry  classes,  it  nevertheless  is  useful  to  clearly  define  the  coordinate  system  conventions 
and  rigorously  develop  the  transfonnations  employed  at  ARL  when  we  provide  instrumentation 
and  data  collection/processing  support  for  projectile  flight  experiments.  It  is  hoped  that  this 
elementary  and  somewhat  tedious  exposition  will  benefit  interpreters  of  body-fixed  sensor  data 
and  other  readers  by  familiarizing  them  with  the  issues  involved  and  thereby  avoid  the  confusion 
that  we  sometimes  have  observed  with  regard  to  such  data. 
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To  this  end,  three  Cartesian  systems  commonly  employed  in  projectile  test  and  evaluation  are 
defined  and  the  mathematics  necessary  to  transform  among  them  are  given  in  this  report. 
Examples  of  idealized  projectile-fixed  sensor  output  from  projectile  trajectory  simulations  and 
sensor  output  recorded  during  flight  experiments  will  illustrate  this  process.  An  additional 
coordinate  system  extensively  employed  in  computer  modeling  and  simulation  of  projectile 
flight  is  treated  in  Appendix  A. 


2.  Coordinate  Systems 


The  first  coordinate  system  is  convenient  for  describing  the  motion  of  a  gun-  or  tube-launched 
projectile  along  its  trajectory.  This  system  is  right-handed  Cartesian  with  its  origin  at  the  launch 
site.  This  is  referred  to  as  the  “earth-fixed”  system  herein,  and  parameters  in  this  system  are 
subscripted  with  an  “E”.  The  axes  are  defined  as  follows  (see  Figure  1): 

•  The  XE  and  YE  axes  define  a  plane  tangent  to  the  earth’s  surface  at  the  origin. 

•  The  ZE  axis  is  perpendicular  to  the  earth’s  surface  with  positive  downward,  i.e.,  in  the 
direction  of  gravity. 


•  The  XE  axis  is  chosen  so  that  the  centerline  of  the  launcher  is  in  the  XE-ZE  plane. 
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Z  E  is  into  the  page.  Y  E  is  out  of  the  page. 

Figure  1 .  Earth-fixed  coordinate  system. 
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Down-range  travel  is  then  measured  along  the  Xe  axis,  deflection  is  measured  along  the  Ye  axis 
(positive  to  the  right  when  one  is  looking  down  range),  and  altitude  is  measured  along  the  ZE  axis 
(positive  downward).  Designation  of  height  above  the  ground  as  a  negative  Ze  component  of 
projectile  location  may  seem  unnatural  to  some,  but  this  coordinate  system  was  employed  by  a 
former  division  chief  at  ARL,  Dr.  Charles  Murphy,  in  his  seminal  work  on  free-flight  motion  of 
symmetric  projectiles  (Murphy,  1963),  and  this  became  the  norm  within  ARL. 

The  second  system  is  convenient  for  aeroballistic  computations  of  rigid  projectiles’  flights  and 
for  describing  the  locations  and  orientations  of  such  projectiles’  components.  This  system  is 
right-handed  Cartesian  with  its  origin  at  the  center  of  gravity  (c.g.)  of  the  flight  body  (see 
Figure  2).  This  is  referred  to  as  the  “projectile-fixed”  system  herein,  and  parameters  in  this 
system  are  subscripted  with  a  “P”. 


Most  of  the  projectiles  instrumented  by  ARL  are  either  spin-stabilized  or  are  designed  to  have 
some  fin  cant-induced  roll  rate  to  mitigate  the  effects  of  asymmetries.  For  such  flight  bodies,  the 
Xp  axis  of  the  projectile-fixed  coordinate  system  usually  lies  along  the  projectile  axis  of 
symmetry,  i.e.,  the  spin  axis  (with  positive  in  the  direction  of  travel  at  launch).  The  Yp  and  ZP 
axes  are  then  oriented  so  as  to  complete  the  right-handed  orthogonal  system.  Because  the  Yp  and 
Zp  axes  are  collinear  with  projectile  radii  in  the  cross  section  through  the  c.g.,  sensors  whose  axes 
are  parallel  with  this  cross  section  are  commonly  termed  “radial”  sensors.  Similarly,  sensors 
whose  axes  are  parallel  with  the  Xp  axis  are  commonly  termed  “axial”  sensors. 

The  third  coordinate  system  is  commonly  employed  to  specify  locations  on  or  near  the  earth’s 
surface  (i.e.,  north,  east,  and  down).  This  is  referred  to  as  the  “navigation”  system  herein,  and 
parameters  in  this  system  are  subscripted  with  an  “N”.  Thus,  north  =  Xn,  east  =  Yn,  and 
down  =  Zn. 
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3.  Transformations  With  Euler  Angles 


The  specification  of  a  vector  in  any  Cartesian  system  can  be  transformed  into  its  corresponding 
specification  in  any  other  Cartesian  system  in  two  steps.  First,  the  origins  must  be  made 
coincident  through  a  translation.  The  second  step  is  to  perfonn  three  successive  rotations  in  a 
specified  sequence.  The  three  successive  rotation  angles  are  called  Euler  angles.  This 
relationship  is  illustrated  for  the  earth-fixed  and  projectile-fixed  systems  in  Figure  3  where  the 
gray  ball’s  position  in  the  earth- fixed  system  is  given  by  PE  and  its  position  in  the  projectile- 
fixed  system  by  Pp .  If  the  location  as  measured  in  the  first  system  of  second  system’s  origin  at 
the  projectile  c.g.  is  0E  ,  then  the  vector  from  that  origin  to  the  ball,  as  measured  in  the  first 
system,  is  LE  ,  which  can  be  obtained  from  the  translation  equation  PE  =  0E  +  LE  .  We  can  then 
derive  Pp  from  LE  by  perfonning  the  Euler  rotations  necessary  to  obtain  the  orientations  of  the 
Xp,  Yp,  and  Zp  axes. 


Figure  3.  Position  of  a  ball  as  measured  in  two  coordinate  systems. 

After  the  translational  offset  to  collocate  the  origins  of  the  earth- fixed  and  projectile-fixed 
systems  has  been  subtracted,  transformations  of  vectors  from  the  earth-fixed  system  to  the 
projectile-fixed  system  via  Euler  rotations  proceed  as  follows: 

1.  Rotate  the  XE,YE,ZE  system  about  the  ZE  axis  through  an  angle  8,  where  8,  is 
measured  in  the  direction  of  rotation  of  the  head  of  a  right-hand  screw  advancing  along  the  ZE 
axis  (see  Figure  4).  Call  the  resulting  axes  XpY^Z,  in  which  Zt  =  ZE  .  is  the  angle 
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between  the  XE  axis  and  the  projection  of  the  Xp  axis  into  the  XE  |  YE  plane.  Thus,  the 
projectile’s  spin  axis  is  in  the  X;  |  Z,  plane  with  its  nose  in  the  +  X,  \Zl  half-plane. 


Figure  4.  First  Euler  rotation. 


If  the  coordinates  of  a  vector  V  are  denoted  by  the  column  matrix 

in  the  Xl,Y1,Z]  system,  then  the 

UJ 

coordinates  are  related  by  the  equations 


"O 

bE 

in  the  earth-fixed  system  and  denoted  by 

( 

b, 

VCE  ) 

lciy 

aj  =  aE  cos  (s,)+bE  sin  (s.) 
b,  =  -aE  sin(s,)+  bE  cos(s,) 


or,  in  matrix  form, 

V,  =E1Ve 


in  which 


E, 


/  cos(£j )  sin(£j )  0^| 
-  sin(£j )  cos(£j)  0 
V  0  0  1, 


(1) 

(2) 

(3) 


2.  Rotate  the  Xl,Yl,Zl  system  about  the  Yj  axis  through  an  angle  S2  measured  in  the 
direction  of  rotation  of  the  head  of  a  right-hand  screw  advancing  along  the  Yj  axis  (see 
Figure  5).  Call  the  resulting  axes  X2 ,  Y2 ,  Z2  in  which  Y2  =  Y, .  82  is  the  angle  between  the 
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projectile  spin  axis  and  the  Xj  axis  so  that  the  projectile  spin  axis  lies  along  the  X2  axis.  In 
other  words,  X2  =  Xp. 


Figure  5.  Second  Euler  rotation. 


If  the  coordinates  of  V  are  denoted  by 


in  the  X2,  Y2,Z2  system,  then  the  coordinates  are 


related  by  the  equations 


or,  in  matrix  form, 
in  which 


a2  =  a{  cos(S2)-Cj  sin(S2) 
b2  =  bj 

c2  =  al  sin(s2)+Cj  cos(sJ 


V2  =E2V!  =E2  EjVe 


"cos (S2)  0  -  sin (ZZ2 )' 


E,  =  0  1 


s i n  (£*2 )  0  cos(<6\ ) 


(4) 

(5) 

(6) 


3.  Rotate  the  X2,Y2,Z2  system  about  the  X2  axis  through  an  angle  83  measured  in  the 
direction  of  rotation  of  the  head  of  a  right-hand  screw  advancing  along  the  X2  axis  (see 
Figure  6).  Call  the  resulting  axes  X3,Y3,Z3  in  which  X3  =  X2  =  Xp  .  S3  is  chosen  so  that 
Y3  =  Yp  and  consequently,  Z3  =  Zp  . 
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If  the  coordinates  of  V  are  denoted  by 
related  by  the  equations 


A 


VC3  J 


in  the  X3,  Y3 , Z3  system,  then  the  coordinates  are 


or,  in  matrix  form, 
in  which 


a3  -  a2 

b3  =  b2  cos(s3)+ c2  sin(S3) 
c3  =  -b2  sin(83)  +  c2  cos(S3) 


V3=E3V2=E3E2E1Ve 


e3 


"0  0  1  " 
0  cos(S3)  sin(S3) 
v0  -sin(83)  cos(83)y 


(V) 

(8) 

(9) 


Recognizing  that  the  X3 ,  Y3 ,  Z3  system  is  the  same  as  the  projectile-fixed  Xp ,  Yp ,  Zp  system 
and  forming  the  product  matrix  TP  =  E3  E2  E, ,  Equation  (8)  becomes 

VP=TEPVE  (10) 


Further,  the  choices  of  Euler  angles  were  made  in  order  to  correspond  with  angles  commonly 
used  to  describe  projectile  orientation  throughout  a  trajectory.  The  first  Euler  angle,  Sj ,  is  the 

azimuth  component  of  the  projectile  heading  relative  to  the  initial  shot  line.  In  ballistic 
tenninology,  this  is  commonly  called  the  yaw  angle  and  is  denoted  by  \\t .  The  second  Euler 

angle,  S2 ,  is  the  elevation  component  of  the  projectile  heading  or,  in  ballisticians’  parlance,  the 
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pitch  angle,  9  .  The  third  Euler  angle,  8^ ,  is  the  rotation  component  of  the  projectile  orientation 
or  the  roll  angle,  (j>  ,l 

When  we  perfonn  the  matrix  multiplication  and  make  the  substitutions  for  the  Euler  angles,  T,1’ 
is  seen  to  be 

%  T2  V 

tep=  t4  t5  t6  (11) 

VT7  ^8  T9  y 

in  which 

T,  =  cos(0)cos(\|/) 

T2  =  cos(0)  sin(\j/) 

T3  =  -sin(9) 

T4  =  sin(9)  sin(cj))  cos(\(/)  -  cos(()))  sin(\)/) 

T5  =sin(9)sin((j))sin(\(/)  +  cos((j))cos(\(/)  (12) 

T6  =  cos(0)sin((j)) 

T7  =  sin(0)  cos((j))  cos(\)/)  +  cos(cj))  cos(\(/) 

T8  =  sin(0)  cos((j))  sin(\)/)  -  sin(cj))  cos(\)/) 

T9  =  cos(0)  cos((j)) 

The  inverse  transfonnation  of  a  vector  from  projectile-fixed  coordinates  to  earth-fixed 
coordinates  is  given  by 

Ve=(tef)TVp  (13) 

Equations  (11)  and  (13)  can  be  used  to  transform  any  vector  between  the  earth- fixed  and 
projectile-fixed  systems.  Acceleration,  velocity,  thrust,  and  line  of  sight  are  examples  of  vectors 
whose  components  typically  are  needed  in  both  systems.  Before  proceeding  to  examples  of 
transformations  of  this  type  of  sensor  data,  another  critical  characteristic  of  the  relationship 
between  the  earth-fixed  and  projectile-fixed  systems  needs  to  be  considered,  i.e.,  the  projectile- 
fixed  system  is  moving  with  respect  to  the  earth-fixed  system  in  both  translation  and  orientation. 
The  angular  rates  describing  these  orientation  changes  need  to  be  formulated  in  each  system  as 
well. 


'  It  is  important  to  note  that  even  though  the  use  of  \|/,  9,  and  <|>  to  denote  yaw,  pitch,  and  roll  is  common  in  the 
literature,  not  every  author  defines  each  of  these  angles  as  is  done  here.  Since  there  is  no  great  agreement,  care  must 
be  taken  always  by  authors  to  rigorously  define  the  systems  they  employ  and  by  readers  to  understand  the  authors’ 
systems. 
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4.  Angular  Rates  in  Projectile  and  Earth  Systems 

The  angular  velocity  of  the  projectile-fixed  system  with  respect  to  the  earth- fixed  system  is 
described  in  the  projectile-fixed  system  by  a  vector  Qp  =  (coXp  coYp  coZp  )T  in  which  coXp  is  the 

angular  velocity  of  the  Yp  |  Zp  axes  about  the  Xp  axis;  coYp  is  the  angular  velocity  of  the 
Zp  |  Xp  axes  about  the  Yp  axis;  and  coZp  is  the  angular  velocity  of  the  Xp  |  Yp  axes  about  the 
Zp  axis.  Each  of  these  velocities  is  defined  herein  to  be  right-hand  positive.  In  ballisticians’ 
tenninology,  this  vector  is  denoted  by  QP  =  (p,q,r)T . 

The  angular  velocity  vector  of  the  projectile-fixed  system  with  respect  to  the  earth- fixed  system 
is  also  given  by  the  sum  of  the  derivatives  of  the  Euler  angles  8, ,  S2 ,  and  S3 .  However, 
recalling  that  each  of  these  Euler  angles  was  defined  in  a  different  coordinate  system,  their 
derivatives  need  to  be  transformed  into  the  projectile-fixed  system  before  corresponding 
components  can  be  added.  With  the  ballistic  terminology  for  the  Euler  angles  adopted  in  the  last 

section,  8,  =(0,0,\j/)T  in  the  X^Y^Zj  coordinate  system,  S2  =(o,0,o)T  in  the  X2,Y2,Z2 

coordinate  system,  and  S3  =  (<j),0,o)T  in  the  X3,Y3,Z3  coordinate  system.  Thus,  QP  can  also  be 
written  in  terms  of  the  Euler  angles  as 

(0)  (0) 

QP  =E,E0  0  +E,  9  +  0  (14) 

3  2UJ  3(oJ  (oj 

Making  the  substitutions  for  E2  and  E3  from  Equations  (6)  and  (9)  and  doing  the  algebra  yields 

^<j)-\j/sin(0)  N 

QP  =  0cos(cj))  +  \j/cos(0)sm((j))  .  (15) 

-  0  sin((j))  +  \j/  cos(0)  cos(cj))^ 

The  relationships  between  the  body- fixed  components  of  projectile  angular  motion  and  the  roll, 
pitch,  and  yaw  rates  are  thus  given  by 

p  =  <j>  -  \j/  sin(0) 

q  =  0  cos(cj))  +  \j/  cos(0)  sin(cj))  (16) 

r  =  -0  sin(4>)  +  \j/  cos(0)  cos(c))) 

Solving  for  the  Euler  angle  derivatives  in  terms  of  the  projectile-fixed  rates  gives 
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4>  =  p  +  [  q  sin(4>)  +  r  cos(<J))  ]  tan(0) 
0  =  q  cos(()))  -  r  sin(4>) 

\j/  =  (q  sin(cj))  +  r  cos(<j)))/ cos(0) 


(17) 


As  seen  later  in  this  report,  there  are  commonly  used,  body-fixed  sensor  systems  designed  to 
measure  p,  q,  and  r  and  other  systems  that  respond  to  <j>,  0,  and  \j/ .  Failure  to  appreciate  the 
differences  among  these  quantities  can  lead  to  erroneous  interpretations  of  sensor  data. 
Conversely,  appreciation  of  these  differences  sometimes  can  provide  additional  insights 
concerning  projectile  dynamics. 


5.  Body-Fixed  Sensor  Locations  and  Orientations 


The  location  of  a  sensor  on  a  flight  body  can  greatly  affect  its  output.  In  some  cases,  the  quantity 
being  measured  is  location  dependent,  e.g.,  acceleration  at  a  point  on  the  body.  In  other  cases, 
the  sensor  response  is  affected  by  something  besides  the  primary  stimulus  that  it  is  designed  to 
measure,  e.g.,  the  gain  of  an  angular  rate  sensor  might  be  g-sensitive.  In  both  instances,  the 
locations  and  the  orientations  of  the  sensors  must  be  known  in  order  to  correctly  interpret  their 
output  in  the  projectile-fixed  coordinate  system  even  before  a  transformation  into  the  earth- fixed 
navigation  system.  Seldom,  if  ever,  will  it  be  possible  to  install  sensors  at  arbitrary  locations  on 
projectiles.  Although  greater  possibilities  might  exist  in  a  developmental  or  test  configuration, 
tactical  configurations  usually  have  only  limited  space  at  predetennined  locations  available  for 
on-board  sensors  and  electronics.  Linear  accelerometers  provide  an  illustrative  example  of  the 
effects  of  sensor  location  and  orientation. 


The  basic  equation  relating  position  (P),  velocity  (V),  and  acceleration  (A)  is 

t  t 

P(t)  =  P(0)  +  V(0)t  +  J*J*A5t.  For  projectile  locations  along  a  trajectory  as  measured  in  the 

0  0 

earth-fixed  system,  the  vector  version  of  this  relationship  is 


'*(t)El 

f  x(0)e^ 

f x(0)E t" 
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fax 

ae 
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(18) 


Beginning  with  acceleration  and  velocity  components  in  the  projectile-fixed  coordinate  system, 
this  relationship  becomes 
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(19) 
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In  ballisticians’  tenninology  the  velocity  vector  in  the  projectile-fixed  coordinate  system  is 

.  Thus,  ballisticians  rewrite  Equation  (19)  as 


f  u) 

V 

= 

r  xpa 
yP 

and  the  acceleration  vector  is 

V 

vwV 

lzp ) 

(wy 
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(20) 


Considering  only  the  axial  acceleration  component,  u ,  the  process  required  to  obtain  this 
quantity  from  projectile-fixed  sensor  data  can  be  complicated.  Designating  a  location  on  board  a 
rigid  projectile  with  unchanging  mass  properties  as  (LXp  ,LYp  ,LZp  )T  where  each  component  is 

measured  from  the  projectile’s  c.g.,  the  ballisticians’  equation  for  the  axial  acceleration  at  that 
point  is  given  by 

AXp  =  u  +  qw  -  rv  -  (q2  +  r2  )LXp  +(pq-r)LYp  +(pr  +  q)LZp  (21) 

Because  gravity  cannot  be  sensed  by  an  accelerometer  on  board  a  projectile  in  free  flight,  the 
idealized  output  (sXp )  of  an  axial  accelerometer  at  this  location  is  thus 

SXp  =u  +  qw-rv  +  (-q2 -r2)LXp +(pq-r)LYp +(pr  +  q)LZp +gsin(0)  (22) 


In  order  to  simplify  notation,  sensed  acceleration  components  parallel  to  the  projectile-fixed  axes 
are  usually  referred  to  as  Acc_I  (+XP  direction),  Acc_J  (+YP  direction),  and  Acc_K  (+ZP 
direction).  Similarly,  sensor  offsets  are  usually  given  as  Ai,  Aj,  and  Ak.  Thus,  Equation  (22)  is 

usually  written  as 

Acc_  I  =  u  +  qw  -  rv  +  (-  q2  -  r 2  )Ai  +  (pq  -  r)  Aj  +  (pr  +  q) Ak  +  g sin(0)  (23) 

The  difference  between  axial  acceleration  (u)  and  an  idealized  axial  accelerometer’s  output 
(Acc_I)  can  be  better  appreciated  graphically.  With  a  six-degree-of-freedom  computer  model, 
the  trajectory  of  an  M483  artillery  projectile  was  simulated  for  a  launch  at  25°  quadrant  elevation 
and  340-m/s  muzzle  velocity.  Figure  7  shows  the  axial  acceleration  throughout  the  trajectory 
and  Figure  8  shows  the  output  of  an  ideal  axial  accelerometer  at  the  projectile  c.g.,  L=(0,0,0)T  . 
The  difference  between  the  two  curves  is  primarily  attributable  to  the  absence  of  the  gravity 
component  in  the  sensed  acceleration  in  Figure  8  because  |qw  -  rv|  <  0.0003  g . 
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Time  Since  Muzzle  Exit  (s) 

Figure  7.  Axial  acceleration  (u)  for  modeled  M483  trajectory. 


Figure  8.  Sensed  axial  acceleration  at  the  c.g.  for  modeled  M483 
trajectory. 


When  the  axial  accelerometer  is  installed  in  the  fuze,  as  is  typical  in  flight  experiments,  the 
location  would  be  approximately  0.45  m  forward  of  the  c.g.  The  sensed  axial  acceleration  at  this 
location,  L  =  (0.45,0,0)T  ,  is  shown  in  Figure  9.  For  an  accelerometer  2  mm  off  the  projectile 
spin  axis,  e.g.,  L  =  (0.45,0. 0014, 0.0014)1 ,  the  sensed  acceleration  is  shown  in  Figure  10. 
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Time  Since  Muzzle  Exit  (s) 

Figure  9.  Sensed  axial  acceleration  0.45  m  forward  of  the  c.g. 
for  modeled  M483  trajectory. 
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Figure  10.  Sensed  axial  acceleration  0.45  m  forward  of  the  c.g. 
and  2  mm  off  the  spin  axis  for  modeled  M483 
trajectory. 

Despite  the  lack  of  resemblance  of  the  curve  in  Figure  10  to  the  curve  in  Figure  7,  axial 
acceleration  (it)  can  be  derived  from  accelerometer  output,  given  initial  conditions,  knowledge 
of  the  location  of  an  axial  accelerometer,  and  knowledge  of  the  angular  velocity,  (p,  q,  r)T  ,  and 
acceleration,  (p,q,r)T  ,  of  the  projectile.  These  rates  are  estimated  by  other  on-board  sensors. 

In  the  preceding  examples,  the  accelerometers  were  assumed  to  be  exactly  situated  at  known 
points  and  to  have  their  sensitive  axes  oriented  exactly  parallel  to  the  projectile  spin  axis. 
Neither  assumption  is  realizable  in  an  actual  system  because  of  location  and  alignment 
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uncertainties  arising  from  tolerances  in  the  manufacturing  and  assembly  processes.  Additionally, 
linear  accelerometers  typically  have  some  level  of  sensitivity  to  accelerations  in  directions  other 
than  the  direction  of  the  measurement  axis.  These  uncertainties  must  be  resolved  in  order  to 
correctly  interpret  the  sensor  output  and  to  derive  desired  quantities  from  the  output.  We  can 
accomplish  this  resolution  by  calibrating  the  individual  sensors  either  in  the  laboratory  or  in  flight. 


6.  Sensor  Locations  in  the  Aeroballistic  Diagnostic  Fuze 


The  aeroballistic  diagnostic  fuze  (DFuze)  (U.S.  Patent  6349652)  was  developed  at  ARL  as  a 
high-g,  projectile-borne  system  for  non-intrusive  measurements  of  projectile  diagnostics  and 
aerodynamic  performance.  The  current  nominal  DFuze  instrumentation  configuration 
(designated  as  the  DF2K1  series)  combines  patented  optical  sensors,  a  unique  constellation  of 
accelerometers  and  magnetic  sensors,  analog  and  digital  electronics,  and  telemetry  components. 
DFuzes  can  and  have  been  built  to  conform  to  individual  projectile  form  factor  requirements. 
The  DF2K1  DFuze  maintains  the  interface  characteristics  of  a  standard  North  Atlantic  Treaty 
Organization-compatible  155-mtn  artillery  fuze.  An  example  is  shown  in  the  upper  half  of 
Figure  11.  Two  of  the  four  solar  light-indicating  transducers  (SLITs)  (U.S.  Patent  5909275) 
used  in  the  solarsonde  system  for  measuring  projectile  heading  and  projectile  crossing  rate  with 
respect  to  the  sun  can  be  seen  in  this  view.  As  the  schematic  in  the  lower  half  of  Figure  1 1 
shows,  the  radial  axes  of  the  projectile-fixed  coordinate  system  are  defined  so  that  SLIT  No.  1  is 
in  the  +XP|Zp  plane.  Locations  within  the  DFuze  body  of  all  the  components  are  specified, 
and  laboratory  calibrations  are  used  to  better  measure  installed  locations  and  orientations 
of  the  sensing  elements.  The  locations  are  then  transformed  to  the  projectile-fixed  system  by 
the  addition  of  the  components  of  the  offset  of  the  DFuze  from  the  projectile  c.g.  For  rigid 
projectiles  with  unchanging  mass  properties,  this  is  a  fixed  offset.  For  flexible  projectiles, 
projectiles  with  components  in  relative  motion  (e.g.,  de-spun  nose),  and  projectiles  with 
changing  mass  properties  (e.g.,  through  fuel  consumption  or  cargo  ejection),  c.g.  can  change  in 
flight,  and  location-sensitive,  projectile-fixed  sensors’  output  will  be  affected. 

Besides  the  four  solarsonde  SLITs,  a  DF2K1  DFuze  contains  six  accelerometers,  a  temperature 
sensor,  and  a  magnetometer  mounted  on  an  ARL-designed  board.  This  board  is  indexed  to  the 
DFuze  housing  so  that  the  sensor  locations  are  duplicated  in  each  DFuze  within  manufacturing 
and  installation  tolerances.  Figure  12  shows  the  forward  side  of  the  board  as  installed  within  the 
DFuze  with  the  projectile-fixed  axes  and  identifying  letters  superimposed  on  the  picture.  The 
component  labeled  A  is  a  3 -axis  magnetometer  with  its  axes  nominally  oriented  parallel  to  the 
projectile  axes  and  with  the  same  parity  as  the  projectile  axes.  The  output  of  some  magne¬ 
tometers  (and  other  non-inertial  sensors)  is  affected  by  the  inertial  forces  to  which  the 
magnetometers  are  subjected  while  they  measure  the  primary  stimulus.  Because  of  this,  the 
location  parameters  as  well  as  the  orientation  parameters  of  these  sensors  must  be  known  in 
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order  to  correctly  interpret  their  output.  In  order  to  simplify  notation,  each  magnetometer’s  three 
axes  are  usually  referred  to  as  Mag_I  (+XP  direction),  Mag_J  (+YP  direction),  and  Mag_K  (+ZP 
direction).  The  component  labeled  “B”  is  a  single-axis  accelerometer  used  to  measure  thrust  and 
drag.  This  sensor’s  axis  is  ideally  situated  on  the  projectile’s  XP  axis  and  oriented  in  the  +XP 
direction  and  is  commonly  called  Acc_I.  The  axial  component  (Ai)  of  this  sensor’s  location  is 

its  offset  from  the  c.g.  of  the  assembled  projectile  with  fuze.  Component  “C”  is  a  temperature 
sensor,  which  can  be  used  to  compensate  for  the  variations  in  bias  and/or  scale  factor  with 
temperature,  which  are  characteristic  of  many  sensors.  The  remaining  unlabeled  components  are 
various  supporting  electronics. 


+XP 

SLIT  #4 

+ Y  P  Axis  is 
into  the  slide. 

Figure  11.  DF2K1  DFuze  configuration. 


Figure  12.  DF2K1  sensor  board  viewed 
from  forward  looking  aft. 
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On  the  aft  side  of  the  board  (see  Figure  13),  there  are  five  accelerometers,  labeled  D,  E,  F,  G,  and 
H.  Accelerometer  D  is  a  bi-axial  device  whose  sensitive  axes  are  ideally  situated  on  the 
projectile  spin  axis  and  oriented  parallel  to  the  Yp  and  Zp  axes.  The  axial  location  of 
accelerometer  D  is  measured  from  the  c.g.  of  the  projectile  with  fuze  assembly.  These  sensor 
axes  are  usually  called  Acc_J  (-Yp  direction)  and  Acc_K  (+Zp  direction).  The  orientation  of  the 
Acc_J  axis  in  the  -Yp  direction  is  an  artifact  of  accelerometer  D  being  mounted  on  the  aft  side  of 
the  sensor  board.  A  parity  correction  of  the  Acc_J  output  is  made  during  data  reduction  to 
provide  acceleration  measurements  in  the  projectile-fixed  coordinate  system.  Accelerometers  E, 
F,  G,  and  H  are  single-axis  devices  whose  output  is  combined  on  board  the  projectile  to  provide 
an  estimate  of  the  projectile  spin  rate.  This  collection  of  sensors  is  called  the  accelerometer  spin 
ring  or  the  accel  ring  for  short.  In  ballisticians’  notation,  the  equations  for  sensed  acceleration  at 
a  point  in  the  +Yp  and  the  +Zp  directions  are 

SYp  =  v  +  ru-pw  +  (pq  +  f)LXp  +(-p2-r2)LYp  +(qr-p)LZp  -  gcos(0)sin(^))  (24) 

SzP  =  w  +  pv  -  qu  +  (pr  -  q)LXp  +  (qr  +  p)L Yp  +  (-  p2  -  q 2  )LZp  -  g cos(0) cos(<|>)  (25) 

With  the  simplified  notation  of  Equation  (22),  these  equations  become 

Acc_J  =  v  +  ru-pw  +  (pq  +  f)Ai  +  (-p2  -r2)Aj  +  (qr-p)Ak-geos(0)sin((t))  (26) 

Acc_K  =  w  +  pv-  qu  +  (pr  -  q)Ai  +  (qr  +  p)Aj  +  (-p2  -  q2  )Ak  -  g cos(0) cos(c)))  (27) 


AFT  Side  of  Board 


+ZP 


+XP  Axis  is  into  the  slide. 


Figure  13.  DF2K1  sensor  board  viewed  from  rear  looking 
forward. 
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Before  any  flight  experiment,  Equations  (22),  (24),  and  (25)  or  alternatively,  (23),  (26),  and  (27) 
need  to  be  evaluated  for  anticipated  projectile  dynamics  in  order  to  determine  required 
accelerometer  characteristics  such  as  dynamic  range  and  frequency  response.  The  natures  of  the 
accelerations  to  which  the  projectile-fixed  sensors  will  be  subjected  and  the  natures  of  the 
acceleration  components  in  an  earth-fixed  or  navigation  coordinate  system  can  be  very  different. 
The  sensor  characteristics  necessary  to  obtain  projectile-fixed  measurements  of  sufficient 
accuracies  to  meet  required  earth-fixed  measurement  accuracies  are  the  design  performance 
requirements.  These  requirements  cannot  be  understood  without  an  appreciation  of  the 
relationship  between  a  sensor’s  location  and  orientation  and  that  sensor’s  stimuli. 


7.  Estimating  Projectile  Spin  Rate 


Spin  rate  is  an  important  characteristic  of  any  projectile’s  flight.  This  rate  (p  in  ballisticians’ 
notation)  is  defined  in  the  projectile-fixed  system  as  the  angular  velocity  of  the  Yp-  Zp  axes 
about  the  XP  axis.  The  DFuze  does  not  include  any  sensor  that  directly  measures  spin  rate. 
Because  of  size,  high-g  survivability,  and  cost  requirements  and  the  necessity  for  sometimes 
estimating  spin  rates  of  300  Hz  and  more,  traditional  angular  rate  measurement  devices  such  as 
gyroscopes  have  not  been  included  in  the  DFuze.  Perhaps  there  will  someday  be  micro-electro- 
mechanical  devices  suitable  for  inclusion  in  future  DFuzes,  but  at  present,  spin  rate  estimates  are 
made  from  accel  ring,  solarsonde,  and  magnetometer  measurements. 


7,1  Accel  Ring  Measurements 

The  four  accelerometers  in  the  accel  ring  (see  Figure  13)  are  installed  in  opposed  pairs  about  the 
projectile  spin  axis  with  the  sensor  axes  of  each  pair  co-linear  with  the  diameter  upon  which  they 
are  situated  and  co-directional  in  orientation.  That  is,  the  sensitive  axes  of  accelerometers  E,  F, 
G,  and  H  are  oriented  in  the  +YP,  +YP,  +ZP,  and  +ZP  directions,  respectively.  The  differences 
between  the  sensed  accelerations  from  two  such  pairs  on  the  YP  axis  (F,E)  and  the  ZP  axis  (H,G) 
are 


Acc_  JF  -  Acc_  JE  =  (AjF  -  AjE)(-p2  -r2) 
Acc_Kh  -Acc_Kg  =(AkH  -AkG)(-p2  -q2) 


(28) 


If  AjE  =  |AjF|  =  AkG  =  AkH  =  R,  then 


Acc_  JF  -Acc_  JE  +  Acc_Kh  -Acc_Kg  =  4R 


'  2  ,  r2+q2^ 


p-  +- 


Acc  _  JF  -  Acc  _  JE  +  Acc  _  Kh  -  Acc  _  KG 

4R 


f  .  r2  +  q2  ^ 
V  1  J 


(29) 

(30) 
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This  arithmetic  is  done  on  board  the  DFuze,  and  the  telemetered  measurement  is  the  quantity 

2  2 

r  +  q 

given  by  Equation  (30).  When  p2  »  — — — ,  as  is  often  the  case,  a  good  estimate  of  the 
projectile  spin  rate  is  obtained  from  the  accel  ring  output  by 

P* 


2  ,  2  \ 

2  r  +  q 

P"+— r*- 


v 


(31) 


7.2  Solarsonde  Measurements 

The  four  SLIT  optical  sensors  installed  in  the  DF2K1  DFuze  are  designed  to  produce  a 
significant  output  when  aligned  with  the  sun  and  almost  no  output  when  not  aligned.  On  board  a 
spinning  projectile,  these  sensors  provide  a  pulse  train,  which  when  combined  with  calibration 
data,  yields  a  measurement  of  the  angle  between  the  spin  axis  and  the  vector  to  the  sun  (gs, 
called  the  solar  aspect  angle)  and  a  solar  roll  history.  Figure  14  shows  a  sample  of  telemetered 
solarsonde  data  from  a  DFuze  mounted  in  a  stationary  spin  fixture.  The  signals  from  all  four 
SLITs  are  combined  on  board  the  projectile  into  a  single  output  stream.  The  spikes  on  the  graph 
of  this  combined  wavefonn  occur  when  the  respective  SLITs  are  aligned  with  the  sun  during  a 
projectile  rotation.  The  combination  of  gains  and  parities  of  the  four  SLITs  are  varied  so  that  the 
output  from  each  sensor  can  be  identified  (as  is  indicated  by  the  numerals  on  the  plot).  Also 
shown  are  two  simple  period  measurements  that  can  provide  a  solar  roll  rate  estimate. 


400  . 
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Time  (s) 

Figure  14.  Sample  of  Solarsonde  data 
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If  a  projectile’s  spin  rate  (p)  or  solar  aspect  angle  (as)  changes  within  a  roll  cycle,  error  is 
introduced  into  the  period-based  solar  roll  rate  estimate.  In  this  instance,  a  more  accurate  roll 
rate  estimate  is  obtained  by  the  correlation  of  the  solar  aspect  angle  at  each  SLIT  aligmnent  with 
the  roll  angle  at  that  alignment  measured  during  laboratory  calibration.  Summation  of  the 
amount  of  roll  required  for  each  succeeding  SLIT  alignment  yields  an  accumulated  roll  angle 
versus  time  history.  First  and  second  derivatives  give  estimates  of  solar  roll  rate  and  roll 
acceleration,  respectively.  These  estimates  are  more  accurate  than  those  obtained  from  period 
measurements  because  of  the  greater  frequency  of  roll  orientation  measurements.  This 
methodology  is  typically  employed  in  the  processing  of  flight  data. 

7.3  Magnetometer  Measurements 

Figure  15  shows  the  data  from  one  of  the  radial  magnetometer  axes  (Mag  J)  during  the  same 
time  period.  Period  measurements  between  signal  extrema  are  used  to  estimate  projectile  roll 
rate  with  respect  to  the  magnetic  field.  In  contrast  to  the  SLITs  in  the  solarsonde  system,  which 
are  designed  to  be  nearly  impulsive  in  their  responses,  vector  magnetometers  provide  continuous 
response  where  the  amplitude  of  the  output  signal  varies  directly  with  the  projection  of  the  local 
magnetic  field  onto  the  sensor  axis.  With  knowledge  of  the  local  magnetic  field  and  a  sensor  that 
is  accurately  located,  aligned,  and  calibrated,  we  can  directly  derive  projectile  heading  with 
respect  to  the  magnetic  field  from  an  amplitude  measurement.  When,  as  is  usually  the  case  with 
actual  implementations,  any  of  these  assumptions  is  violated,  advanced  processing  techniques 
are  employed  to  remove  the  effects  of  these  errors. 


14.00  14.10  14.20  14.30  14.40  14.50 

Time  (s) 

Figure  15.  Sample  of  radial  magnetometer  data. 
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For  these  solar  and  magnetic  sensor  systems,  what  is  being  measured  is  not  the  projectile  spin 
rate  but  the  rate  at  which  the  optical  and  magnetic  sensors  are  swept  past  the  sun  and  the  earth’s 
magnetic  field,  respectively.  The  equations  for  these  rates  are  the  same  as  the  equation  for  the 
derivative  of  the  Eulerian  roll  angle  (Equation  [17])  with  appropriate  change  of  variables.  With 
S  and  M  subscripts  to  indicate  solar  and  magnetic,  respectively,  these  equations  are,  with  the 
necessary  changes, 

<j)s  =  P  +  [qsin((j)s)  +  rcos(cj)s)  ]tan(0s)  (32) 

<j>M  =  P  +  [  q sin(4>M )  +  r  cos(cj)M)  ]tan(0M)  (33) 


Since  solar  roll  rate  and  magnetic  roll  rate  are  entirely  analogous,  further  discussion  is  restricted 
to  solar  roll  rate.  With  some  manipulation,  Equation  (32)  becomes 


4>s  =  P  +  [qsin(^s)  +  rcos(c^)] 

=  qsin(4>s)  +  r  cos(c[)s) 
COS(0s) 

=  p  +  [\j/seos(0s)]  tan(0s) 


COs(^>  I  tan(0s) 
cos(es)  1 


cos(0s)  tan(0s) 


(34) 


With  the  recognition  that  the  term  [\j/s  cos(0s)]  is  the  component  of  projectile  yawing  rate 

perpendicular  to  the  plane  defined  by  the  projectile  spin  axis  and  the  solar  vector,  that 
0S  =  90  -  as,  and  the  realization  that  Equation  (34)  holds  for  any  earth-fixed  field,  the  relationship 
between  projectile  roll  rate  with  respect  to  any  such  field  and  projectile  spin  rate  can  be  stated 
thus.  Roll  rate  is  the  sum  of  spin  rate  and  the  product  of  the  out-of-plane  yawing  rate  and  the 
tangent  of  the  complement  of  the  angle  between  the  spin  axis  and  the  field  vector.  Roll  rate  and 
spin  rate  are  equal  when  there  is  no  yawing  motion  (\j/  =  0)  or  the  spin  axis  is  perpendicular  to 
the  field  vector  (a  =  90, .'.  0  =  0)  and  unequal  when  \j/  ^  0  and  a  ^  90  .  This  relationship  and  the 
relationship  between  spin  rate  and  accel  ring  data  have  sometimes  provided  crucial  insight  into 
projectile  dynamics  not  attainable  from  other  measurements. 


8.  Summary 


Body-fixed  sensors  combined  with  a  telemetry  system  provide  measurements  of  projectile  in¬ 
flight  dynamics  not  achievable  with  stationary  or  remote  sensing  systems  such  as  radar  or 
cameras.  Body-fixed  sensor  systems  on  board  the  projectiles  provide  continuous  data  throughout 
a  flight,  but  many  of  the  desired  performance  metrics  are  defined  in  an  earth-fixed  coordinate 
system.  It  is  therefore  crucial  that  the  relationship  between  these  two  systems  be  understood  for 
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the  proper  design  of  projectile-fixed  sensor  systems  and  for  the  correct  interpretation  of  sensor 
data  obtained  during  flight  experiments. 

The  mathematics  describing  this  relationship  have  been  presented  herein  and  examples  have 
been  provided  for  the  sensor  systems  included  within  the  DFuze  measurement  system,  i.e., 
accelerometers,  vector  magnetometers,  and  solarsondes.  Because  these  mathematics  hold  for 
any  vector  definable  in  these  coordinate  systems,  this  work  is  applicable  to  any  sensor  system 
measuring  vector  quantities,  e.g.,  angular  rate  sensors,  velocimeters,  gimbaled  seekers,  etc. 


21 


9.  Reference 


1.  Murphy,  C.H.  Free  Flight  Motion  of  Symmetric  Missiles;  BRL-1216;  Ballistics  Research 
Laboratory:  Aberdeen  Proving  Ground,  MD,  July  1963. 


22 


Appendix  A.  Plane-Fixed  Coordinate  System 


The  basic  force  and  moment  equations  from  which  the  free  flight  equations  of  motion  of  a  rigid 
projectile  are  derived  are  simplified  when  formulated  in  a  non-rolling  coordinate  system  with  the 
same  X  axis  as  the  projectile-fixed  system.  The  plane-fixed  coordinate  system  is  such  a  non¬ 
rolling  system  with  the  added  stipulation  that  the  plane-fixed  Y  axis  is  in  the  horizontal  plane. 
Thus,  the  plane-fixed  Z  axis  is  in  the  vertical  plane  containing  the  spin  (X)  axis.  Designating  the 
axes  with  an  FP  subscript,  the  plane-fixed  system  is  shown  in  Figure  A-l. 


Y 


X  FP  into  the  page. 


Figure  A-l.  Plane-fixed  coordinate  system. 


The  plane-fixed  system  is  used  extensively  in  computer  trajectory  simulations  since  the 
differential  equations  of  motion  can  be  evaluated  at  larger  integration  time  steps  in  this 
coordinate  system  than  in  the  projectile-fixed  system.  The  plane-fixed  system  already  was  seen 
in  Section  3  of  this  report  although  it  was  not  so  named  at  that  time.  The  X2 ,  Y2 ,  Z2  system  that 

resulted  from  the  second  Euler  rotation  is  the  plane-fixed  system.  The  transformation  matrix 
from  earth-  to  plane-fixed  vectors,  given  by  substituting  the  Euler  angles  and  evaluating  the 
matrix  multiplication  in  Equation  (5),  is 


^  cos(9)  cos(\|/)  cos(0)  sin(\|/)  -  sin(9)^ 
-  sin(\|/)  cos(\|/)  0 

vsin(0)cos(\|/)  sin(0)sin(\|/)  cos(0) 


(A-l) 


The  transfonnation  matrix  from  the  plane-fixed  to  the  projectile-fixed  system  has  already  been 
given  in  Equation  (9).  Substituting  the  Euler  roll  angle,  this  becomes 
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1 


'0  0 

TFPP  =  0  cos(()))  sin(cf))  (A-2) 

v0  -  sin(4>)  cos(())) 
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